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Abstract 

We consider an elementary discrete process which starts from purely 
random configuration and leads to well-ordered and stable state. Com- 
plete analytical solution to this problem is presented. 



1 Statement of the problem 

Simple models of evolution may sometimes lead to complicated and unexpected 
behavior, cf. e.g. [6]. We consider the following elementary problem. Let us 
choose a finite set of n points {pi} distributed arbitrarily (in particular ran- 
domly) on the plane. Let us further label these points in a unique but otherwise 
quite arbitrary way: 

Po, pi,..., Pn-l (1) 

Of course, given any set one has n! different ways of labelling its elements. We 
can further treat these points as two-dimensional sequence of vectors attached 
to the origin of the Cartesian coordinates. 

Define next the following "cyclic" sequence of vector differences and denote 
them Pi(l): 

po(l) = pi-po (2) 
= P2-P1 

Pn-2(1) = Pn-l-Pn-2 
Pn-l{l) = P0-Pn-1 

It represents a sequence of n new vectors which, after appropriate shifting, may 
also be attached to the origin of coordinates. We can consider these points as 



the next phase of some discrete time evolution and repeat the process itera- 
tively many times. Note that this process is strictly deterministic, that is, once 
chosen the initial points and their labelling, the rest of the discrete evolution is 
fully determined regardless of how many steps has been performed. The only 
chaotic element in this process is included in the initial distribution, unless this 
distribution is deliberately chosen as regular. 

Each transition t — > t + 1 in this discrete evolution can be written using more 
suggestive matrix notation as: 
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for the x coordinate and similarly for the y coordinate, where t — 0,1,2,... 
Consequently, starting from the initial configuration, after t evolution steps, we 
get: 
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where the power is understood as matrix power, of course. (Matrix which 
appears in ([3]) and ([4]) is a special case of so-called circulant matrix, cf. [5], for 
pedagogical review see [4].) 

Using different, perhaps less evocative notation, this algorithm can also be 
written in coordinate representation as: 
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Equations ([5]) coincide with with the higher order differences [5] . 

In the final section complete analytical solution to this problem will be given 
using Discrete Fourier Transform (DFT), see e.g. [2]. 



2 Numerical experiments 

This problem may easily be investigated numerically. The experiments starting 
from many different random initial distributions give surprising and decidedly 
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counter-intuitive results. Naively, one would expect that this algorithm, when 
starting from purely random distributions, would just produce more random 
distributions. However, it turns out that the final configuration depends simply 
on whether n is even or odd. In case of n even, after roughly n steps of discrete 
evolution, we get exactly two different separate loops of points. One of these 
loops contains points with odd labels whereas the other contains points with 
even labels. On the other hand, in case of n odd we get a single loop of points. 
All configurations are asymptotically point symmetrical with respect to the 
origin. (In other words, the center of mass of the system is in the center of 
coordinates.) 

These results may be qualitatively understood. The crucial thing here is 
that the number of points n is finite and they follow a cyclical pattern. Since 
the t-th order forward derivatives ([5]) involve t+1 terms, when t exceeds n some 
numbers Xi(0) appear more than once. When t 3> n then each Xi(0) appears 
many times in a linear yet complicated combination. 

All figures below were made using Mathematica (cf. Fig. 5). We also made 
some animated gifs which illustrate this phenomenon and are available from the 
authors. 

3 Analytical solution 

The statement of the problem was given by the first author. He also performed 
the numerical experiments described above. The complete analytical solution 
which we shall present below was given by the second author. 

We may consider the coordinates separately. The calculations below concern 
the x coordinate but may be applied to y (as well as to other possible coordinates 
in higher dimensions). 

Let us denote initial values of coordinates by 

x o {0),xi(0), . . . ,x n -i(0). 

Recall that the discrete time evolution is described by the equations 

Xj(t + 1) = £Ej_|_i(t) — Xj(t), t — 0, 1, ... ,n — 2 
x n -x{t + 1) = x {t) - x n -i(t). 

The evolution is linear, therefore it may be diagonalized using discrete Fourier 
transform 

n— 1 

x k {t) = - r Y, x i® u ~ ik > fc = 0,l,...,n-l, (7) 

V n j=Q 

where oj — exp(27ri/n). Fourier transform of the eq. ([6]) reads 

x k {t+l) = (uj k -l)x fe (i). (8) 

The general formula is 

x k (t) = (Lu k -l)'i fe (0) (9) 
3 
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The inverse Fourier transform of eq. (JHJ) is 

1 n — 1 ^ n — 1 n— 1 

u(t) = T ri l (t)/ = -V( w i -i)yT Ij (o) U - ji (io) 

v k=0 fc=0 j=0 

1 n— 1 n— 1 

= -^ f {u k -l)*u lk ^ f x j (0p- ih / = 0,l,...,n-l. 

n fc=0 j=0 

From this point we will write equations for two coordinates. 

In the case of odd n, absolute value of \ui k — 1| reaches maximum of r = 
2cos(7r/2n) for k = (n ± l)/2, u k — — exp(±7ri/n). Hence for large t we have 
the following approximation 

xAt) w r'(-l)' +t cos — — .A + sin — — B (11) 

\ n n J 

y \ t , 7r(Z + t/2) . 7r(Z + i/ 2 ) ^ 
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where A,B,C,D are numerical coefficients: 

A = - y7-l) J 'cos — 2^(0) (12) 
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We note that the following relations hold 

xi{t + 2) = ~r 2 x l+1 {t) (13) 
Vi{t + 2) = -r 2 y l+1 (t). 

We have in this case 

(:;SH ( - i,,+, (c B D )(zt) <»' 

where 

= 

n 

Hence for large t and AD — BC ^ 0, the points Po(t),Pi{t), ■ ■ ■ ,p n -i(t) lie on 
an ellipse 

(C 2 + D 2 )x 2 - 2{AC + BD)xy + (A 2 + B 2 )y 2 = {AD - BC) 2 r 2t . (15) 
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In the case of even n, absolute value of \uj k — 1| reaches maximum value equal 
2 for k = n/2, ui k = — 1 / Hence for large t we have an approximate relation 
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More detailed investigations, when higher eingenvalue is taken into account, 
give the following asymptotics: 
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where r\ = 2cos(7r/n) and A\, Bi, C\ 1 D\ are defined similarly to (fT2"]) with 27r 
instead of 7r. 



4 Conclusion 

The process presented in this paper may easily be generalized to more spatial 
dimensions. The method described in the previous section applies also in this 
case. Just as in the 2D case the final configuration depends solely on whether 
n is odd or even (cf. Fig. 4). 

Let us finally note that this example of discrete time evolution leading from 
random to ordered distribution complies with several informal requirements 
formulated by some distinguished mathematicians when they spoke about aes- 
thetics or even beauty of their discipline. Namely, given result should be easily 
stated, yet counter-intuitive, deeply non-obvious, surprising. It should have rig- 
orous and elegant proof. Finally - as G. H. Hardy used to emphasize - it should 
be perfectly useless [3]. 
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Mathematica code used for generating animated gifs illustrating the 
effect described in this paper 

(* Mathematica code *) 

(* M random labelled points on the plane *) 
M = 25; 

ro := Random [Real, {-1,1}, 10] 
Table [{ro, ro , ro}, {m, 1, M}] ; 
points = 7.; 

(* differences between consecutive *) 
(* vectors generate M new vectors *) 

F[P_] := Table [ 

If[k < Length [P], P[[K+1]]-P[[k]] ,P [ [1] ] -P [ [Length[P] ] ] ] , 
{k, 1, Length [P]} 

] 

(* iteration and normalization *) 

Nest[F, points, i] 

G[i_] := 

Norm [Nest [F, points, i] ] 

(* color of points: odd - green, even - red *) 

b[i_] := If[ 0ddQ[i], 0.4, 0.0] 

Mathemaatica code (continued) 
Odrawing points *) 

H[k_] := Module [ 

{to}, 

Show[ 
Table [ 
Graphics3D [{ 
Hue [b [k+i] ] , 
PointSize [r] , 
Point [ 

{t0[[i]] [[1]] , t0[[i]] [[2]] , t0[[i]] [[3]]}] 
}], 

{i, 1, Length [t0]} 



G 



Background -> RGBColor [0 . 2 , 0.2, 0.3], 
AspectRatio -> 1, 
Axces -> False, 
Boxed -> True, 

PlotRange -> {{-R, +R>, {-R, +R>, {-R, +R», 
Viewpoint -> {100, -50, 50}, 

PlotLabel -> StyleForm ["times" <> ToString[k], FontSize -> 12] 

] 

(* size of dots *) 
r = 0.008; 

(* size of picture *) 
R = 0.13; 

tO = SessionTime [] ; 

Table [H [V [k] ] , {k, 1, 510, 1} ] 

SessionTime [] - tO 

Export ["animation.gif " , 7.7., "GIF", ImageResolution -> 100] 
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Figure 1: Typical discrete time evolution for n = 24 (left) and n = 25 (right). 
For even n two different loops emerge whereas for odd n single loop comes out. 
Since these loops grow with time, here and in all the remaining figures, coordi- 
nates Xi are normalized, i.e. divided by yY^^o ( x i) 2 ( an d similarly for y). 
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Figure 2: 20 cases of initially randomly distributed n = 50 points after 300 
steps of discrete time evolution. Points pk with odd and even labels k (k = 
0,1, 2, ...,7i— 1) are marked with different colors in order to show that in the 
final configuration they belong to different loops and do not mix. At this stage 
there are various symmetric shapes but they all tend to two ellipses, see the 
next figure. 



9 





V 




1 




/} 

v 




1 


• 


- 





/ 


* 


- 





/ 




9 




§ 


- 







9 




m 


\ 


\ 








\ 




m 




o 




m 


• 

• 


o 



Figure 3: The same as in the previous figure but after 1500 steps of discrete 
time. 



10 



Figure 4: 20 cases of initially randomly distributed n —51 points after 300 steps 
of discrete time evolution. Single loop is visible with various degree of entan- 
glement. Note that points with even and odd values of the labelling parameter 
also do not mix. There are many different shapes but they all tend to some 
ellipse, cf. eq. (fT5|) . see also the next figure. 
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Figure 5: The same as in the previous figure but after 1500 steps of discrete 
time. 
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Figure 6: 3D case reveals the same qualitative behavior, n — 24 (left), 
(right). 
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